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ABSTRACT 


A model is suggested in which some of the important features of the 
circulating , 'flow inside the two-dimensional near wake are derived by 
assuming a slow viscous flow. The theory considers the flow away from 
the body base. It is found that there is a region of constaivt speed 
merging, as we go downstream, into a region of stagnation-apex flow. 

The velocity returning from the rear stagnation point along the center 
streamline is shown to be a slowly varying function of the "wedge-angle" 
of the wake and to be roughly one half the velocity at the edge of the 
shear layers driving the wake-cavity flow. These results seem to be in 
agreement with experimental data. 
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Orr.-'t NGft 47-lOZ-OOl while the author was in residence at ICASE, NASA 
Lorgiey Research Center, Hampton, VA 23665. 








1. Introduction 


A complete description of the structure of the two dimensional near wake 
behind a body in motion depends on knowledge of the whole flow field. In 
particular one has to know, or derive, the way the boundary layer separates 
at the base and becomes a shear layer. This shear layer is driven by the 
inviscid free stream and in turn drives the recirculating flo>-< confined inside 
the near- wake cavity. 

Experiments at low to moderate (Re<200) Reynolds numbers [1] indicate 
that the recirculating flow field in the cavity is composed of three parts: 

The flow very near the base, the central region, and the flow field near the 
rear stagnation point where the cavity closes on itself and where the flow, 
driven by the shear layers, turns around towards the base. The data show 
that in the central region the velocity on the stagnation streamline (connecting 
the I ear stagnation point to the center of the base) is constant, 
independent of the Reynolds number and is about 1/5 the free stream value. 

flumer ical studies for supersonic wake and Re=500 {2] indicate? that in the 
near- wake the ’-'iscous and pressure terms dominate while a short distance down- 
stream of the -.-ear stagnation point the two merged shear layers have boundary- 
layer- like St -i!;. t’lre. 

The abov'e observations suggest that one may be able to model the flow in 
the santral and rear stagnation regions of the recirculating flow field by 
considering a slow viscous fluid bounded by converging streamlines which 
"represent" the bounding shear layers. Away from the apex, i.e. the rear 
stagnation point, the velocities along the outer and center streamlines are to 
be constant. Nearer the apex the velocities will depend linearly on the 
distance from the stagnation point and very close to it the variation will be 
parabolic - as suggested by Heimenz (See for example Batchlor [3]). 


In section 2 the Stokes' equation will be solved for the central region, 
and it v;ill be found that the returning velocity alcng the center stream line 
cannot be assigned arbitrarily but that It is dictated by the constraint that 
on the center line the pressure gradient normal to it must vanish. 

In section 3 the linear stagnation region is considered, if is fouiid 
that in it the pressure gradient must be radial. Theref«jre, the solution 
in the constant velocity region found in section 2 and the .olution for this 
stagnation region nan be matched without changing the wedge angle, 

In section 4 v/e examine the final parabolic approacit to stagnation. It 
is found that, as in the central (constant velocity) region. there is a con- 
straint on the returning velocity along the centerline. This condition, 
however matches the one found in section 2 only for limiting small apex 
angles. Thus in general this region will produce a discontinui cy in the 
slopes of the streamlines. However, because our model problem possesses no 
natural reference lengtn, this ^inai stagnation region can be made as small 
2 ' desired and vn i ) not affec.t the solution elsewhere. An interesting result 
'.n this imm-diaL^* aehj:ioorhO‘<(i of the stagnation apex 1S that the boundary 
■ t -eaml MI- 3 S cannot Porn, an apex angln of 26^ = n radians f;LCi i;.e the returning 
seiocitv a’eng the t .nterline (6 - 0; is proportional to sec d. — i.e. it 
as radicb y</tL. 

2. Uie c oi.stant ( radial ) velocity region 

In this reg’on the geometry is as follows: 



FIGURE 1 


The boundary conditions are “"*^0 streamlines w - 0, 0 = 26^ 

and u, = 0 on the same streamlines. Since we assumed the inertia forces 
6 

to be nugatory the equation to be solved is Stokes', which is bi-harmonic 
for the streamfunction, 

= 0 . ( 2 . 1 ) 

With the above boundary conditions, however, we shall get a sink- type flow. 
We must stipulate that along the center .ine the fluid is moving away from 
the apex (which is not in ou»* region of consideration now). Therefore we 
consider the following problem: 


V. 



FIGURE 2 


•t is cler' that m polar coordinates the solution to (2.1) with its attendant 
bo; . !a / conditions must oe of the form 

III = rf(e) 


from which 


\ " r Ifl " 
‘'e " ■ 
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and equation (2.1) for the stream function becomes 


r^y4^ = f(iv) ^ 2f« + f = 0. 


( 2 . 2 ) 


The boundary conditions for the half-wedge (see figure 2) are 


f(0) = 0 


f(0) = -U 


f(6„) . 0 


f(0„) 


(2.3) 


The solution of (2.2) subject to (2.3) is 


f(0) = A^(sin0-0ros0) + B^0sin6 - U^0cos0 


(2.4) 


where 


A. = 


K- !(!!^)9.sin0 V ! 


° e^sinV ru^'^o-'-o^o 


(2.5) 


and 

“ ; 2 “° 2 :" ! 

0 -Sin 6 I 0 ) 

0 0 

If we wish the center line e = 9 q to be a true line of syrmietry for the 
complete wedge it is not sufficient to have Ug(0^) = 
also require that the pressure gradient normal to that streamline vanish; 
i.e. (Sp/30)„ Q = 0. Now 

0-0Q 


iB. = li(f‘'+f ) 
80 r^ ' 


(2.7) 
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whe^e p is the coefficient of viscosity. The symmetry requirecnent leads 
therefore to 

(A^+Uo)sine^ + B^ccse^ = 0 (2.8) 

Since A and B depend on 0„, U„ and V^/U„ we find that for a given 

driving velocity and a given semi-apex angle 6^ equation (2.8) yields 
a value for V^/U^. The variation of = V^/U^ with is shown in 
figure 3 . Note that for 050^11/2, 3 ^ varies monotonically between 0.5 

and 2 /ir = . 637 . Thus for most values of 0 ^ that one expects to encounter, 
in experiments for example, the value of = V^/U^ will be confined to 
fairly narrow bounds. This seems to be in agreement with experiments [ 1 ] 
where at different (low) Reynolds number the same constant value of u^/U^ 

was measured along the stagnation streamline. The measured value of 
!u^/U^l = 0.2 would give our results the right order of magnitude if U^, 
the velocity characteristic of the "bottom" edge of the driving shear-layer, 
is taken to be about 1/3 the local free stream value. 


3 . The (linear) stagnation region 

We assume that from a certain distance, r^, the velocity along the 
boui.Jing streamlines begins to decelerate linearly in r. This is net an 
unreasonable assumption and is in line with Heimenz' results [ 3 ] for stagnation 
point flow. Thus we set in this region (r < r^) 

ij; = r^f(6). 


Consequently 

u^ = rf'(0). 




Ug = -2rf(0), 





and the bi-harmonic equation reduces to 



f(0) 0 f(e^) = 0 

f‘{0) = -u^ f‘(e^) = Vj . 


(3.2) 


Since r has not been scaled yet we shall say that this stagnation region 

is defined by r £ r^ = 1. Then in order for the velocity along the bounding 

streamline to be continuous in the two regions r ^ 1 and n £ 1 we must 

set U, = U and V, = V . 

1 0 1 0 

The solution to the system (3.1), (3.2) is 


f = A^{sin20-20) + B^(cos29-l) - U^8 

where 

U (l+6jtane -20„ 

A = _o . 1 0 0 

1 4 e^-tana„ 

0 0 


(3.3) 


(3.4) 


Uq (sin0pCOse^)(l+3j+20ptan0^)-U+8^)6Q 

“ T ’ ~ " (0 -tan9 )(sin0^cos0„) 

0 0 0 0 


Vi 

Q = 1. — 0. 

% Uj 


(3.6) 


It is easily established that for the present case 3p/30 = 0 everywhere 
and hence there is no constraint on the value of 3|. It follows that it was 
legitimate to require continuity of the radial velocity on the bounding stream- 
lines at r = 1 ; i.e., to set, as we did, 6^^ = 6^. 






2 

It is interesting to note that the choice of » r f(0) is the only one 
leading to 3p/30 = 0 everywhere. This is seen by examining the expression 
for the pressure gradient: 


_L ie = 

pr 39 



^2 30 



If Ug and are calculated from = r‘’^f{9) we find 


(3.7) 


r'"'^(2-m)(r+m2f) (3.8) 

It is seen from (3.8) that 9p/38 vanishes everywhere for m = 2. This 

is fortunate since otherwise we would have f 6^ and would not have been 

able to match the velocities on the center streamline of the regions at 

r = r =1, or at any other station. We are now is a position to combine 
0 

the results of this section with the previous section to specify a flow-field 
of a slow viscous recirculating fluid in a wadge-like region (with the possibl 
exception of the very near neighborhood of the apex): 


r>l 


j 

^ rfA (sin0-0cos0) + B 0sin0 - U„0cos6] 
0 0 0 

1 

O<0<8„ 

,l; = 

1 

1 

(3.9) 


^ r^[Aj(sin28-29) Bj(cos20-1) - U^6] 

0£r£l 


O<0<0„ 

with Aq, and Bj given by (2.5), (2.6), (3.4) and 

(3.5) respectively 


There remains the question of the continuity of and its derivatives, at 
r = 1, at all 0 ■ and not just along the bounding streamlines, 0=0 and 
6 = 0Q, where the boundary conditions wore chosen to assure us of continuity. 
Clearly there is no analytic matching uniformly valid in 6 . Numerical 


calculations of f and f as functions of e for various 6^ are shown 
in figures 4 and 5. It is seen that the matching is quite good. 


4. Final approach to the stagnation ape x 

For O^rf^r^^l we might have a parabolic appraoch to stagnancy, i.e., 


il, = r^f(0) 


(4.1) 


This small r 2 is likely to be of the order of the diffusion length proper 
to this problem, ~ o(y/pU^), where we assumed tacitly that « 1 and 
hence U U L . 

1»R_ = . 

From (4.1) we have 


u = r f ' (e) 

r 


Ug = -3r‘^f(6) 


and the bi -harmonic equation for this case becomes 


rVijj = + 10f'= + 9f = 0 


(4.2) 


with boundary conditions 

f(0) = 0 f(e^) = 0 

f*(0) = -U 2 f'(0o) = 


(4.3) 


where continuity of velocity along the bounding streamline (0=0) demands 
that or = V’^2' (^-2) and (4.3) for f(0) we find 

■j '■^2 

f = A2(sin0 --|sin39) + B2(cos0-cos30) - ^ sin30 


(4.4) 


where 




2 

n e^+pgcose^-i 

stn^0^ 

0 


B 


_^2 . 
2 “ T 


2cos0^-B2 

sinBo 


B 


2 



(4.5) 


(4.6) 

(4.7) 


From (3.3) it follov/s that the symnietry requirement (3p/30 = 0 at 0 = 0^) 
takes the form 


f"(0^) + 9f{0g) = 0 (4.8) 

Using (4.4), (4.5) and (4.6) we find that S 2 is constrained to take the 
value 


^2 "" 2COS0 

0 


(4.9) 


Evidently ^ ’’ range of interest, O<0^^/2 , 

^2^®o^ -- equality taking place only at 6^ = 0. Thus in the 

present model there is no continuous matching possible between the linear and 
parabolic stagnation regions. Either, if we match the velocities, there is 
a discontinuity in the slopes of the bounding and other streamlines, or -- 
if we leave the wedge shape undisturbed — there is a discontinuity in the 
velocity field. This was to be expected in a relatively crude model as the 
one employed here. However, the extent of the parabolic region is very 
small and its effect on the rest of the flow field should be small. The main 
point is that the two wdjor regions — that of constant velocity and the linear 
stagnation one — do match well to yield a flow field with properties charac- 
teristic of available low Reynolds numbers data. 


^ 9 - „ . ..... 


It is interesting to note that the expression for effectively 

excludes the possibility of a blunted wedge shaped cavity with closure by 
streamlines perpendicular to the center line, i.e., 0^ “ ^^2 ’ "fact, 

one would guess that at the apex the closure is cusped. 
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Figure 5.- Matching of the radial velocity distribution -ry- versus 
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